REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  data  sources, 

gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 

of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Service,  Directorate  for  Information  Operations  end  Reports 

1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget 

Paperwork  Reduction  Project  (0704-0188)  Washington,  DC  20503.  , 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY) 

03/03/2014 

2.  REPORT  TYPE 

Interim  Research  Performance  Report 
Quartelry 

3.  DATES  COVERED  (From  -  To) 

10/01/2013-12/31/2013 

4.  TITLE  AND  SUBTITLE 

Adaptive  Meshing  of  Ship  Air-Wake  Flowfields 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

N00014-14-1-0010 


5C.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Dr.  Steve  L.  Karman 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
University  of  Tenneesee  at  Chattanooga 
615  McCallie  Ave 
Chattanooga,  TN  37403-2504 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Same 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


12.  DISTRIBUTION  AVAILABILITY  STATEMENT 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT  ” 

The  objective  of  the  proposed  research  is  to  deliver  enhanced  mesh  adaptation  capabilities  that  account  for  the  chaotic 
unsteady  nature  of  the  flowfield  in  the  landing  approach  path.  The  research  will  explore  three  distinct  meshing  methods 
to  handle  the  dynamic  aspect  of  this  case.  The  three  methods  include  a  hierarchical  Cartesian  hexahedral  method,  an 
all-tetrahedral  mesh  method  and  a  physics-based  point  placement  meshless  method.  The  initial  phase  of  the  research 
will  develop  pseudo-steady  analysis  procedures  for  each  method,  with  appropriate  mesh  adaptation  capabilities.  The 
goal  of  the  research  is  to  compare  and  contrast  the  adaptive  meshing  technologies  of  the  three  approaches. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

SAR  1 

18.  NUMBER 

OF  PAGES 

6 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Dr  Steve  Karman 

REPORT 

u 

h.  ABSTRACT 

S 

1  c.  THIS  PAGE 

JL 

19b.  TELEPONE  NUMBER  (Include  area  code) 

(423)  425-5492 

Standard  Form  298  (Rev.  8-98) 
Prescribed  by  ANSI-Std  Z39-18 


INSTRUCTIONS  FOR  COMPLETING  SF  298 


1.  REPORT  DATE.  Full  publication  date, 
including  day,  month,  if  available.  Must  cite  at 
lest  the  year  and  be  Year  2000  compliant,  e.g., 
30-06-1 998;  xx-08-1 998;  xx-xx-1 998. 

2.  REPORT  TYPE.  State  the  type  of  report,  such 
as  final,  technical,  interim,  memorandum, 
master's  thesis,  progress,  quarterly,  research, 
special,  group  study,  etc. 

3.  DATES  COVERED.  Indicate  the  time  during 
which  the  work  was  performed  and  the  report 
was  written,  e.g.,  Jun  1997  -  Jun  1998;  1-10  Jun 
1 996;  May  -  Nov  1 998;  Nov  1 998. 

4.  TITLE.  Enter  title  and  subtitle  with  volume 
number  and  part  number,  if  applicable.  On 
classified  documents,  enter  the  title  classification 
in  parentheses 

5a.  CONTRACT  NUMBER.  Enter  all  contract 
numbers  as  they  appear  in  the  report,  e.g. 

F336 1 5-86-C-5 169. 

5b.  GRANT  NUMBER.  Enter  all  grant  numbers 
as  they  appear  in  the  report,  e.g. 

1 F665702D1257. 

5c.  PROGRAM  ELEMENT  NUMBER.  Enter  all 
program  element  numbers  as  they  appear  in  the 
report,  e.g.  AFOSR-82-1234. 

5d.  PROJECT  NUMBER.  Enter  al  project 
numbers  as  they  appear  in  the  report,  e.g. 
1F665702D1257;  ILIR. 

5e.  TASK  NUMBER.  Enter  all  task  numbers  as 
they  appear  in  the  report,  e.g.  05;  RF0330201 ; 
T4112. 

5f.  WORK  UNIT  NUMBER.  Enter  all  work  unit 
numbers  as  they  appear  in  the  report,  e.g.  001 ; 
AFAPL30480105. 

6.  AUTHOR(S).  Enter  name(s)  of  person(s) 
responsible  for  writing  the  report,  performing  the 
research,  or  credited  with  the  content  of  the 
report.  The  form  of  entry  is  the  last  name,  first 
name,  middle  initial,  and  additional  qualifiers 
separated  by  commas,  e  g.  Smi:h,  Richard,  Jr. 

7.  PERFORMING  ORGANIZATION  NAME(S) 
AND  ADDRESS(ES).  Self-explanatory. 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER.  Enter  all  unique  alphanumeric  report 
numbers  assigned  by  the  performing 
organization,  e.g.  BRL-1234;  AFWL-TR-85- 
4017-Vol-21-PT-2. 

9.  SPONSORING/MONITORS  AGENCY 
NAME(S)  AND  ADDRESS(ES).  Enter  the  name 
and  address  of  the  organization(s)  financially 
responsible  for  and  monitoring  the  work. 

10.  SPONSOR/MONITOR'S  ACRONYM(S). 

Enter,  if  available,  e.g.  BRL,  ARDEC,  NADC. 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S).  Enter  report  number  as  assigned 
by  the  sponsoring/  monitoring  agency,  if 
available,  e.g.  BRL-TR-829;  -215. 

12.  DISTRIBUTION/AVAILABILITY 
STATEMENT.  Use  agency-mandated 
availability  statements  to  indicate  the  public 
availability  or  distribution  limitations  of  the  report 
If  additional  limitations/restrictions  or  special 
markings  are  indicated,  follow  agency 
authorization  procedures,  e.g.  RD/FRD, 
PROPIN,  ITAR,  etc.  Include  copyright 
information. 

13.  SUPPLEMENTARY  NOTES.  Enter 
information  not  included  elsewhere  such  as 
prepared  in  cooperation  with;  translation  of; 
report  supersedes;  old  edition  number,  etc. 

14.  ABSTRACT.  A  brief  (approximately  200 
words)  factual  summary  of  the  most  significant 
information. 

15.  SUBJECT  TERMS.  Key  words  or  phrases 
identifying  major  concepts  in  the  report. 

16.  SECURITY  CLASSIFICATION.  Enter 
security  classification  in  accordance  with 
security  classification  regulations,  e  g.  U,  C,  S, 
etc.  If  this  form  contains  classified  information, 
stamp  classification  level  on  the  top  and  bottom 
of  this  page. 

17.  LIMITATION  OF  ABSTRACT.  This  block 
must  be  completed  to  assign  a  distribution 
limitation  to  the  abstract.  Enter  UU  (Unclassified 
Unlimited)  or  SAR  (Same  as  Report).  An  entry  in 
this  block  is  necessary  if  the  abstract  is  to  be 
limited. 


STANDARD  FORM  298  Back  (Rev.  8/98) 


ONR  Grant  Award  No.  N000 14-14-1  -00 1 0 

Adaptive  Meshing  of  Ship  Air-Wake  Flowfields 

Progress  Report 
March  3,  2014 


The  objective  of  the  proposed  research  is  to  deliver  enhanced  mesh 
adaptation  capabilities  that  account  for  the  chaotic  unsteady  nature  of  the 
flowfield  in  the  landing  approach  path.  The  research  will  explore  three 
distinct  meshing  methods  to  handle  the  dynamic  aspect  of  this  case.  The 
three  methods  include  a  hierarchical  Cartesian  hexahedral  method,  an  all- 
tetrahedral  mesh  method  and  a  physics-based  point  placement  meshless 
method.  The  initial  phase  of  the  research  will  develop  pseudo-steady 
analysis  procedures  for  each  method,  with  appropriate  mesh  adaptation 
capabilities.  The  goal  of  the  research  is  to  compare  and  contrast  the  adaptive 
meshing  technologies  of  the  three  approaches. 

Common  configurations,  including  carrier  and  aircraft,  will  be  simulated 
with  each  method.  Where  feasible  common  software  components  will  be 
used  by  all  three  methods.  These  will  include  geometry  libraries,  spacing 
field  libraries,  flux  functions,  time  integration  routines  and  others.  Three 
distinct  computer  codes  are  being  developed  in  the  process  and  will  link 
with  these  common  components.  The  creation  of  the  three  computer  codes 
has  begun.  Initial  results  for  two  of  the  methods  are  provided  in  this  report; 
the  hierarchical  Cartesian  mesher/solver  and  the  all-tetrahedral  solver. 

The  Cartesian  code,  named  OctFlow,  has  capabilities  similar  to  the 
SPLITFLOW  code  developed  by  Lockheed  Martin.  [1]  OctFlow  utilizes  an 
Octree  data  structure  and  performs  cut-cell  operations  at  geometry 
boundaries.  A  second-order  spatial  finite-volume  scheme  has  been 
incorporated  with  explicit  first  order  backward  time  integration.  The 
Cartesian  mesh  is  generated  automatically  based  on  the  input  geometry, 
which  is  supplied  as  a  triangulated  surface  mesh.  The  cells  intersected  by  the 
geometry  are  handled  using  the  “cut-cell”  approach,  which  is  basically 
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creating  arbitrary  polyhedral  elements  with  appropriate  surface  boundary 
conditions.  Any  cells  completely  outside  the  computational  domain  are 
tagged  external  and  not  solved  in  the  flow  solution.  This  preliminary 
capability  does  not  yet  include  solution-based  adaptation. 

The  all-tetrahedral  code,  named  TetFlow,  reads  an  externally  generated 
mesh  comprised  of  tetrahedral  volume  elements  and  triangular  surface 
elements.  Meshes  for  this  code  are  currently  generated  using  Pointwise.[2] 
This  code  also  uses  a  second  order  spatial  finite-volume  scheme  with  first 
order  explicit  time  integration.  Once  the  basic  flow  solution  capability  is 
developed  the  insertion  of  additional  bodies  (i.e.  aircraft)  will  be 
incorporated  into  the  code  using  the  mesh  rupture  technique  in  three 
dimensions.  [3] 

Both  codes  use  the  same  flux  function  to  compute  the  inviscid  flux  on  the 
control  volume  boundaries.  The  current  formulation  is  the  compressible 
inviscid  Euler  equations  in  conservation  form.  The  codes  are  being  modified 
to  use  a  low-speed,  preconditioning  scheme  to  handle  the  mix  of  flow 
regions  anticipated  during  the  research.  [4] 

A  notional  aircraft  carrier  has  been  simulated  with  the  two  codes  and  is 
shown  below.  The  surface  mesh  from  the  3D  mesh  generated  by  Pointwise 
serves  as  the  geometry  for  the  OctFlow  code.  A  crinkle  cut  view  of  the 
meshes  from  the  two  codes  is  shown  in  Figure  1.  The  cell-cutting  process  in 
OctFlow  uses  the  triangulated  surface  mesh  to  cut  into  the  volume  elements. 
As  a  result  a  finer  surface  tessellation  is  generated,  as  can  be  seen  in  the  left 
side  of  the  image.  The  sizes  of  the  Cartesian  volume  elements  are  dictated 
by  the  size  of  the  surface  triangles,  so  the  local  element  sizes  are  comparable 
between  the  meshes.  A  top  view  of  the  surface  meshes  from  each  code  is 
shown  in  Figure  2.  A  magnified  view  of  the  two  meshes,  shown  in  Figure  3, 
reveals  just  how  different  the  two  meshes  are  given  the  same  boundary 
tessellation. 


Some  preliminary  solutions  have  been  computed  on  the  two  meshes.  The 
flow  conditions  are  Mach  0. 1  inviscid  flow  directed  head-on  to  the  carrier. 
Figure  4  shews  a  side-by-side  comparison  of  pressure  contours  on  the  carrier 
surface.  The  contour  levels  are  set  to  match  exactly  between  the  two  meshes. 
The  element  counts  are  quite  different  between  the  two  meshes.  The 
TetFlow  mesh  contains  515,939  nodes  and  2,904,133  elements.  The 
OctFlow  mesh  contains  1,313,044  nodes  and  342,761  elements.  One  factor 
contributing  to  the  larger  node  count  in  the  OctFlow  mesh  is  the  surface 
cutting  process.  And  tetrahedral  meshes  always  have  larger  element  to  node 
ratios  compared  io  other  element  types,  especially  hexahedrai  elements. 


The  solution  runs  times  are  dramatically  different.  These  initial  results  do 
not  include  any  adaptive  refinement  in  the  OctFlow  mesh,  which  would 
increase  the  element  counts.  The  TetFlow  mesh  size  should  remain  constant 
as  the  strategy  involves  inserting  bodies  into  this  mesh  via  the  rupture 
technique.  Adaptation  of  the  TetFlow  mesh  will  be  through  point  movement, 
not  refinement. 

A  next  step  in  the  research  is  to  complete  the  modifications  to  low-speed 
preconditioning  flux  scheme  and  incorporate  in  both  OctFlow  and  TetFlow. 
The  third  code,  named  PointFlow,  is  also  in  development.  It  is  obviously 
more  complicated  than  the  other  two.  Previous  work  in  3D  needs  to  be 
combined  with  knowledge  gained  in  2D  by  graduate  student  Philip  Fackler 
in  a  code  framework  consistent  with  OctFlow  and  TetFlow.  Some 
components  will  be  directly  shared,  such  as  geometry  and  spacing  libraries. 
Others  will  be  dramatically  different,  such  as  post-processing  visualization 
files. 


References 


1.  Karman,  S.  L.,  Domel,  N.  D.,  "SPLITFLOW:  Progress  in  3D  CFD 
with  Cartesian  Omni-Tree  Grids  for  Complex  Geometries",  AIAA- 
2000-1006,2000. 

2.  "Pointwise."  Pointwise,  Inc.,  http://www.pointwise.com. 

3.  O'Connell,  M.,  and  Karman,  S.,  "Mesh  Rupturing:  A  Technique  for 
Geometry  Insertion  and  Significant  Mesh  Movement",  AIAA-2013- 
0148,  51st  Aerospace  Sciences  Meeting  Including  the  New  Horizons 
Forum  and  Aerospace  Exposition,  Grapevine,  TX,  2013. 

4.  Reed,  C.  L.,  Karman  S.  L.,  "Implementation  of  Low  Speed 
Preconditioning  in  the  Splitflow  Code",  AIAA- 1997- 1867,  1997. 


